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ABSTRACT 


The  numerical  evaluation  of  the  theoretical  variance/mean  ratio  of  the  neutron 
population  in  the  subcritical  condition  of  the  AFRRI-TRIGA  reactor  is  given  for  (1)  the 
case  of  no  delayed  neutron  groups;  (2)  the  case  of  one  equivalent  delayed  neutron  group; 
and  (3)  the  case  of  six  delayed  neutron  groups.  The  Courant-Wallace  equations  were 
solved  to  obtain  the  theoretical  variance/mean  ratio  using  the  mean  prompt  neutron 
lifetime  of  39  x  10-®  seconds  in  the  AFRRI-TRIGA  reactor. 
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I.  INTRODUCTION 


The  mean  number  of  neutrons  and  delayed  neutron -forming  precursors  in  a  re¬ 
actor  may  be  estimated  by  deterministic  equations  but  the  actual  number  at  time  t  will 
be  different.  The  instantaneous  number  of  neutrons  and  precursors,  both  possessing 
a  discrete  rather  than  a  continuous  probability  distribution  function,  will  fluctuate  at 
random  about  the  mean.  A  statistical  parameter  which  may  be  used  to  measure  the 
fluctuations  is  the  variance/mean  ratio.  References  2  and  6,  through  the  use  of  the 
probability  balance  and  two  probability  generating  functions,  derive  equations  describ¬ 
ing  the  stochastic  processes  within  the  reactor's  core  in  terms  of  the  first  and  second 
moments  of  the  neutron/precursor  probability  distributions.  These  equations  may  be 
solved  for  the  variance/mean  ratio.  It  is  the  purpose  of  this  report  to  numerically 

evaluate  the  theoretical  variance/mean  ratio  of  the  neutron  population  for  three  cases 
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using  the  fission  data  for  U  from  references  3  and  5  and  the  mean  prompt  neutron 

—3  1 

lifetime  in  the  AFRRI-TRIGA  reactor  which  is  approximately  39  x  10  seconds.  The 
symbolism  of  reference  6  will  be  followed.  Case  1  is  when  it  is  assumed  that  there 
are  no  delayed  neutron  groups.  Case  2  is  when  there  is  one  equivalent  delayed  neutron 
group.  Case  3  is  for  six  delayed  neutron  groups. 
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II.  GLOSSARY  OF  BASIC  SYMBOLS 


k  -  k!  +  S  c-  =  the  mean  number  of  neutrons  formed  (both  instantaneous  and  de- 
i  1 

layed)  per  neutron  lost  in  the  reactor  through  absorption  or  escape. 
k!  =  the  mean  number  of  neutrons  formed  instantaneously  per  neutron  lost  in  the 
reactor  through  absorption  or  escape. 


c  =  L  c{  . 

i 

c-  —  the  mean  number  of  radioactive  neutron-forming  precursors  of  type  i  formed 
per  neutron  lost  in  the  reactor  through  absorption  or  escape. 

0Nn  “  variance  of  N  =  -  N  . 

®Nj  -  NCj  -  N  Cj  . 

0ij  -  CiCj  -  C|  Cj  . 

Cj  =  the  total  number  of  precursors  of  type  i  in  the  reactor. 

€i  = 

Xj  =  decay  constant  of  precursor  type  i. 

T  =  mean  prompt  lifetime  of  a  neutron  in  the  reactor. 

f  =  f(x,y2?y2>  .  .  . )  =  f-type  probability  generating  function 

=  S  S  S'...  xn  y-.ml  y2m2  .  .  .  p(n,m1,m2,  .  .  .)  . 
n  m2  1  z 

p(n,m^Jm2,  .  .  .)  =  the  probability  that  the  absorption  or  escape  of  one  neutron  in  the 

reactor  will  lead  to  the  formation  of  n  instantaneous  neutrons  and  m^  precur¬ 
sors  of  type  i. 

n  =  number  of  neutrons  formed  instantaneously  per  neutron  lost  in  the  reactor 
through  absorption  or  escape. 
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number  of  precursors  of  type  i  formed  per  neutron  lost  in  the  reactor  through 


absorption  or  escape. 


it  is  to  be  noted  that  k=  n  +  L  rm  since  kT  n  and  c-T)ci  “  m  ”  £  m- . 

i  1  1 

F  =  F(x,y1,y2,  .  .  .  ;t)  =  F-type  probability  generating  function 

=  Z  Z  Z  ...  xNyiCly2C2  ...  P(N,C1,C«,...;t). 

N  Cl  C2 

P(N,  Cj_,  C2,  .  .  .  ;t)  =  the  probability  that  N  neutrons,  C-^  precursors  of  type  1,  C2  pre¬ 
cursors  of  type  2,  etc.  ,  are  present  in  the  reactor  at  time  t. 
t  =  time . 

N  =  the  total  number  of  neutrons  in  the  reactor. 

subscripts  on  f  and  F  -  the  subscripts  on  f  and  F  indicate  partial  differentiation,  i.e.  , 


<  )  -  the  angle  brackets  indicate  evaluation  of  the  enclosed  term  at  the  point 

x=yri. 

N  -  the  bar  over  a  symbol  indicates  the  mean  value,  thus,  N  =  the  mean  value 
of  N. 

6jj  =  Kronecker  delta 


6y  =  0  when  i  f  j 


6|j  =  1  when  i  =  j. 


the  probability  of  achieving  a  fission  event  per  neutron  lost  through  absorp¬ 


tion  or  escape. 
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V  -  the  number  of  neutrons  formed  per  fission  event. 

/3  =  the  mean  number  of  delayed  precursors  formed  per  fission  neutron  produced. 

=  the  mean  number  of  delayed  precursors  of  type  i  formed  per  fission  neutron 
produced. 

prime  symbol  -  the  prime  symbol,  except  for  kT,  is  a  shortened  notation  to  indicate 

division  by  N,  i.e.  ,  0,xrNJ  =  . 

1NJN  N 

kg  “  mean  square  of  the  number  of  neutrons  formed  per  neutron  lost  through  ab¬ 
sorption  or  escape. 

=  1-k+c. 

=  half -period  of  radioactive  precursor  of  type  i. 

<x  =  percent  yield  of  p/3  precursors  of  type  i. 
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III.  GENERAL  EQUATIONS 


The  equations  in  references  2  and  6,  describing  the  fluctuations  of  the  neutron/ 


precursor  probability  distributions,  are 


(1  k+c)  0NN  -  .2^  q  0M  (f^)  +  (1-k+c) 


N 


(1) 


which  results  in  one  equation  as  we  let  i  =  1,2, 3, 4,  5,  and  6,  and 

6 

~Ci0NN+  (1-k+C+Cj)  0m  -  Cj  - 
which  results  in  six  equations  for  the  six  values  of  i,  and 


<fix)  -2Cij  N 


(2) 


— ci  0 Nj  _cj0Ni+  +  ei>  0ij  +  2ciJ  N  (3) 

which  results  in  21  equations  when  i  =  1,2, 3, 4,  5,  and  6,  and  when  j  =  1, 2, 3, 4,  5,  and 
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6  since  the  mixed  partial  derivatives  of  0,  namely,  0-  j  =  0i  :  ,  are  equal  when 

A  JA  JB  XB 

iA  =  iB  and  when  jA  =  JB- 


IV.  CASE  1:  NO  DELAYED  NEUTRONS 


The  assumption  of  no  delayed  neutron  groups  means  we  assume  that  c  =  =  0, 

m^  =  0,  and  that  k  -  k?.  Equations  (1),  (2),  and  (3)  reduce  to 

(1_k>0NN  =  <2-fxx  +  l~k)  N  .  (4) 

With  m^  =  0,  the  f-type  probability  generating  function  becomes 

f(x)  =  L  Xn  p(n)  (5) 

n 

and 

f(x)xx  "fxx  C§x)  "  E  n(n-l)x  '  p(n)  (6) 


and 


=  E  n(n-l)  p(n) 
n 

2 

-  E  n  p(n)  -Snp(n) 
n  n 


Using  (10),  equation  (4)  rearranges  to 

^NN  =  1+  n2  -  n 
N  2  (1— k) 


(7) 

(8) 
(9) 

(10) 


(11) 


By  definition,  k,  the  effective  multiplication  factor  of  the  reactor  system,  is  the 
mean  number  of  neutrons  formed  per  neutron  lost  in  the  system  through  escape  or  ab¬ 
sorption.  Thus 


5 


k=iv 


(12) 


where  4  =  the  probability  of  achieving  a  fission  event  per  neutron  lost  through  absorp¬ 
tion  or  escape  and  where  v  =  the  number  of  neutrons  formed  per  fission  event. 

Assuming  that  4  and  v  are  two  independent  sets  (events), 

k  =  iv  ~  i  v  =  iv  (13) 

since  by  definition  4  =  £  • 

Since  k  =  k’ ,  therefore 

k  =  4p  =  n  .  (14) 

Thus 

?  =  (4^)2  =  i2  vZ  .  (15) 

Substitute  (15)  into  (11),  and  obtain 

0NN  =  1  +  42  P2  -  k  .  (16) 

N  2(l-k) 

By  selecting  values  for  the  dependent  variable  k  we  may  obtain  the  variance/mean 
ratio  of  the  neutron  population  for  various  conditions  of  subcriticality  for  the  no  delayed 
neutron  case  using  equation  (16).  Table  II  (page  17)  and  Figures  1  and  2  (page  18) 
summarize  the  numerical  results. 

V.  CASE  2:  ONE  DELAYED  NEUTRON  GROUP 

To  obtain  equations  for  one  equivalent  delayed  neutron  group  we  let  cj  =  c.  Equa¬ 
tions  (1),  (2),  and  (3)  become 

(1-k+c)  0NN  -  e  0NC  =  ( \  <  W  +  1_k+c)  N  (17) 

•  ~c  0NN  +  (1_k+c+e)  $NC  ~  *  0CC  =  (<~fxy^  "  2c>  N  (18> 
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(19) 


-c0NC  +  e0CC  =  (!(fyy)+C)N 

where 

0NC  =  NC  -  hf  C* 

and  where 

0CC=  C2  -  C  . 

Add  (17),  (18),  and  (19)  and  obtain 

(1-k)  (0NN  +  0NC )  =  (‘F^xx'1  +  +  2  ^fyy^  +  1  -  k)  N  . 

Since 

(fxx>  =  £  n(n-l)p(n,m) 

**  n,m 


and 


and 


and 


<fxv>  =  £  nm  p(n,  m) 

n ,  m 


<f  )  =  £  m(m-l)  p(n, m) 
yy  n,  m 


(ko-k)  =  £ 

^  n,m 


2  _| 

(n+m)  -  (n+m)  J  p(n,m) 


(20) 

(21) 


(22) 


(23) 


(24) 


(25) 


(26) 


it  is  seen  that 

2  ^2“^)  =  2  ^xx)  +  +  2^fyy^  ' 

Substitute  (27)  into  (22)  and  obtain 

(1-k)  (0NN  +  0NC^  ~  2  ^2-k)  +  (l~k)  j  N  .  (28) 

Solve  (17),  (19),  and  (28)  for  0NN/Nby  Cramer’s  rule  (see  reference  6  for  de¬ 
tails)  and  obtain 
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(29) 


Ann  =  1  +  0*2 -k)  (i-k+c) 

N  2  (1-k)  (1-k+c+e ) 


— c' 
2  C 


l-k+c+e 


where 


2  f xy  +  f yy 


x=y=l 


=  2  nm  +  m2  -  m  . 


(30) 

(31) 


By  definition,  c  is  the  mean  number  of  precursors  formed  per  neutron  lost 
through  absorption  or  escape.  Thus 

c  =  k/3  (32) 

where  fS  =  the  mean  number  of  delayed  precursors  formed  per  fission  neutron  produced. 

For  the  AFRRI-TRIGA  reactor 

e  =At  (33) 

-  0.077  x  39  x  10~6  (34) 

s:  3 . 0  x  10-6  (35) 

where  0.077  sec-1  is  the  equivalent  one-group  decay  constant  computed  from  the 
6-group  thermal  fission  data  in  reference  5  by  the  following  expression 


^equivalent  6 


(36) 


1 -group 


Ij  \ 

i=l  Ai 


By  definition  m  is  the  number  of  delayed  neutrons  formed  per  neutron  lost.  Thus 
c  =  m  =  Xv  (3  •  (37) 

It  is  assumed  that  X,  v,  and  (3  are  statistically  independent.  Hence 

(38) 


and 


m2  =  (  )2  =  42  /32  =  42  ,2  £2 


(39) 
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By  the  same  logic 


nm  =  (  iv  (1-/3)  x  iv  /3  ) 


(40) 


=  i2  v2  /3  -  £.2  v2  /S2  . 


2  ,  ,2  o2 


(41) 


Substitute  (38),  (39),  and  (41)  into  (31)  and  obtain 
c’  -  2  £2  v2  P  -  i  v  B 

2 

where  the  second-order  term  involving  B  is  neglected. 


(42) 


It  is  realized  that  the  assumption  of  statistical  independence  between  p  and  /3  is 
an  approximation  and  that  the  result  of  this  assumption  is  to  state  in  the  averaging 
process  that  the  mean  of  the  product  of  p  and  /3  is  equal  to  the  product  of  the  separate 
means  when  in  actuality  the  mean  of  p/3  should  involve  sums  of  terms  like  p/3  times  an 
unconditional  probability  and  a  conditional  probability.  Equation  (38)  states,  in  essence, 
that  we  are  considering  £  and  @  as  constants  and  are  combining  the  fluctuations  of  the 
neutron  and  precursor  formations  into  the  p  value.  Overall,  the  value  of  the  last  term 
in  (29),  the  only  term  involving  c*  and  hence  (38),  is,  in  any  case,  small,  and  thus  to 
a  first  order  approximation  the  assumption  of  statistical  independence  between  p  and 
/3  introduces  only  a  negligible  error  which  is  of  a  magnitude  less  than  the  standard 
deviation  of  the  final  ratio  obtained  by  propagating  the  precision  indexes  of  the  experi¬ 
mental  data  for  the  various  other  parameters  in  (29). 

Substitute  (13),  (37),  and  (42)  into  (29)  and  obtain 


g2  p2  (1-2/3)  -  k 
2  (  (1-k)  +  k/S) 


l  + 


(43) 


N 


where  it  is  to  be  noted  in  arriving  at  the  solution  that  the  constant  e  is  of  negligible 
importance . 


v 
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Table  II  (page  17)  and  Figures  1  and  2  (page  18)  summarize  the  0N^/  N  ratios  for 
various  values  of  k. 

VI.  CASE  3:  SIX  DELAYED  NEUTRON  GROUPS 

Twenty -eight  simultaneous  equations  resulting  from  the  expansions  of  (1),  (2),  and 
(3)  must  be  solved  to  obtain  the  0^^/  N  ratio  in  the  case  of  six  delayed  neutron  groups. 
These  28  equations  are 


,0NN~ 

'C10N1"C20N2  "f30N3_c40N4 

e50N5  c6 

52?’- 

II 

I(f  > 
2X  xx; 

+  AJ 

(44) 

;1%N 

+  (A  + 

cl)0Nl  _  e#ll  ~  £2^12  “  c3013  “  c4014  ' 

-C5015' 

~ e  6016 

=  <flx>  ' 

-  2c^ 

(45) 

:2%N 

+  (A  + 

C2)%2  "  f  l^l  “  e2022 

3023~«4024 

5‘ 

■ce^e 

=  <f2x>  " 

■  2c2 

(46) 

:3%N 

+  (A  + 

e3^0N3  "  e  1031  “  c2032  “  c 

3033  ~e4034' 

~c5035' 

"c60 36 

=  <f3x>' 

-  2c3 

(47) 

40NN 

+  (A  + 

c4^ 0N4  “  e  1041  “  e 2042  ”  e 

3043  -e4044' 

"c5045 ' 

"e6046 

=  <f4x>- 

-  2c4 

(48) 

o0NN 

+  (A  + 

O  0N5  “  c 105 1 "  c2052  ~  e 

3053  _c4054' 

_e5055' 

~e6056 

=  <^f5x^  ' 

~2c5 

(49) 

60NN 

+  (A  + 

c 6^0N6  “  e  1061  “  e2062  “  e 

3063 “  e4064 

_e5065' 

"e6066 

=<w 

-  2c  6 

(50) 

"2c10N1 

+  2e10^_1  = 

<fll>  4 

-  2c^ 

(51) 

_C10N2  “  C20N1  +  (e2 

+  ex)  0i2  = 

<f12> 

(52) 

“°1^N3  “  C30N1  +  (e3 

+  ?i) 0i3  = 

(f13> 

(53) 

-C10N4  ~  C40N1  +  (e4 

+  Ci)  014  = 

<*14> 

(54) 

~C10N5  “  C50N1  +  (e5 

+  ex) 015  = 

<f15> 

(55) 

_C10N6  “  C60N1  +  (c6 

+  eiKe  = 

^  f16^ 

(56) 
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-2c20n2  +  2e2?i22  _^f22^  +  2c2  (57) 

_C2^N3  "  c3%2  +  ^3  +  e2  ^23  =  ^23^ 

_C2(^N4  _  C4^N2  +  (f4  +  e2  )^24  =  ^24^  (59) 

_C2^N5  "  C5^N2  +  (f5  +  e2  ^25  =  ^f25^  (60) 

~C2^N6  "  C6^N2  +  (*6+  e2  ^26  =^f26^  (61) 

-2c3^N3  +  2e3^33  =  ^33^  +  2c3  (62) 

-c30N4  “  C4^N3  +  (f4+  e3^34  =  ^f34^ 

~C3^N5  “  c5%3  +  (e5  +  e3  )035  =  ^35^ 

~C3^N6  “  c6%3  +  ^e6  +  c3  ^36  =  ^36^  (65) 

“^C4^4  +  ^^4^44  =  ^44)  +  ^c4  (66) 

-C4I^N5  "  C50N4  +  <c5  +  c4^045  =  ^f45^ 

~C40N6  "  C6^N4  +  <e6  +  e4^46  =  ^f46^  (68^ 

— 2c5%5  +  2e5%5  =  ^f55^  +  2c5  (69) 

“c5%6  "  C6^N5  +  <c6  +  c5>056  =<^f56^  (70) 

-2c60n6  +  2^6^66  ”  ^66^  +  2c6  C71) 

where  A  =  1-k+c  (72) 


and  where  the  prime  mark  over  the  0  indicates  division  by  N,  i.e.  , 

0'  .*m  .  (7: 

NN  n 

235 

The  constants  used  in  solving  (44)  through  (71)  were  obtained  from  '  U  neutron 

induced  fission  data  in  reference  3  and  from  the  thermal  fission  delayed  neutron  data 
235 

for  U  in  reference  5  and  are  listed  below. 
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V  =  2.47  ±  0.03  (±  standard  deviation)  (74) 

V2  =  7. -32  +  0.15  (75) 

=  0.0064  +  0.0003  .  (76) 

Substitute  (32)  into  (72)  and  obtain 

A  =  1-k+c  (77) 

=  1.0  -  0.9936  k  .  (78) 

By  definition 

€.  =X.T  (79) 

l  l 

since 

In  2 

X-  —  rp  (SO) 

1 


where  is  the  half -period  as  described  on  page  473  of  reference  4. 
Substitute  (80)  into  (79)  and  obtain 


_  r  In  2 


ci  =  ^  ,t 


) 


(81) 


The  six  values  of  from  reference  5  are  55 . 72 ,  22.72,  6.22,  2.30,  0.610,  and 
0.230  sec  for  delayed  neutron  groups  i  =  1  through  i  =  6  respectively.  Inserting  these 


values  into  (81)  we  obtain 


«1  = 

0.48515  x  10 "6 

(82) 

^2  = 

1.18982  x  10-6 

—6 

(83) 

e3 

4.34610  x  10 

-6 

(84) 

H  = 

11.75336  x  10 

-6 

(85) 

c5 

44.31597  x  10 

-6 

(86) 

f  6  = 

117.53365  x  10 

(87) 
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By  definition  ci  is  the  mean  number  of  precursors  of  type  i  formed  per  neutron 


lost  through  absorption  or  escape.  Thus  by  analogy  to  (32) 

ci  =  k0i  (88) 

where  is  the  mean  number  of  delayed  precursors  of  type  i  formed  per  fission  neutron 
produced. 

Thus  (88)  may  be  computed  by  the  expression 
kc^  ko^ 

ci  =  2.47  x  100  =  24?  (89) 

where  a ±  is  the  percent  yield  of  v$  of  precursors  of  type  i. 

Substitute  the  six  values  of  from  reference  5,  namely,  0.052,  0.346,  0.310, 
0.624,  0.182,  and  0.066,  into  (89)  and  obtain 


c^  -  2. 112  x  10  ^  x  k  (90) 

C2  =  14. 016  x  10  ^  x  k  (91) 

-4 

Cg  -  12.544  x  10  xk  (92) 

c^  =  25.280  x  10  ^  x  k  (93) 

Cg  =  7.360  x  10  ^xk  (94) 

-4 

Cg  -  2. 688  x  10  x  k  .  (95) 

By  (23)  and  (40) 

<fxx>  =  n(n~!)  =  n2  -  n  (96) 

=  i2u2  (1-0)2  -  tv  (1-/3)  •  (97) 

By  definition  and  by  analogy  to  (40) 

<fix>  =  iiip 

=  [  Zv&i  ]  [  iv  (1-8)  1  (98) 
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(99) 


=  £2  p2  j8i  (1-8)  • 

Similarly  when  i  f  j 

(fij>  =  mjmj  =  (Sl'jSj)  (10°) 

=  (ioi) 

and  when  i  -  j 

(fij)  =  <fji>  =  mi(mi-i)  (i°2) 

=  IWO  C^i3i-1]  <103) 

=  t2v2fS2  -  iv Pi  ■  (104) 


The  assumption  of  statistical  independence  between  p  and  /3  is  used  in  the  evalua¬ 
tion  of  <fxx>,  and  (^j)  • 

Substitute  (13),  (74),  (75),  (76),  and  (89)  into  (97),  (99),  (101),  and  (104)  and 


obtain 

|<fxx)  =  0.59226  k2  -  0.49680  k  (105) 

(flx)  =  k2  x  2.51780  x  10-4  (106) 

<f2x>  =  k2  x  16.70903  x  10_4  (107) 

<f3x)  =  k2  x  14.95420  x  10_4  (108) 

(f^)  =  k2  x  30. 13730  x  10~4  (109) 

<f5x)  =  k2  x  8.77415  xlO'4  (110) 

<f  >  =  k2  x  3.20447  x  10~4  (111) 

(f^)  =  k2  x  10  8  x  5.35185  -  2. 112  x  10  k  (112) 

<fl2>  =  k2  x  10"8  x  35.51682  (113) 

<fl3)  =  k2  x  10"8  x  31. 78675  (114) 
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<f14>  =  k2  x  10-8  x  64.06003  (115) 

<f15)  =  k2  x  10-8  x  18.65039  (116) 

(fir)  =  k2xl0'8x  6.81145  (11?) 

<f22>  =  k2  x  10"8  x  235.70258  -  14.016  x  10~4  k  (118) 

<f23>  =  k2  x  10“8  x  210.94843  (H9) 

(f24)  =  k2  x  10“8  x  425. 12567  (i20) 

<f9r>  =  k2  x  10-8  x  123.77076  (i21) 

(f26)  =  k2  x  10-8  x  45.20324  (122) 

<f__)  =  k2  x  10-8  x  188.79297  -  12.544  x  10~4  k  (123) 

Ou 

(f34)  =  k2  x  10-8  x  380.47765  (i24) 

<f„r>  =  k2  x  10“8  x  110.77197  (125) 

o  o 

(f36)  =  k2  x  10"8  x  40.45585  (i26) 

<f  >  =  k2  x  10~8  x  766.77905  -  25.280  x  10  4  k  (127) 

<f'  =  k2  x  10"8  x  223.23210  (i28) 

x  45 

<f46)  =  k2  x  iO"8  x  81.53093  (i29) 

<f  )=k2xl0'8  x  64.99377  -  7.360  x  IO"4  k  (130) 

(f5g)  =k2xl0_8  x  23.73685  (131) 

(L ,)  =  k2  x  10'8  x  8.66911  -  2.688  x  10“4  k  .  (132) 

The  28  equations,  (44)  through  (71),  were  solved  by  hand  on  a  Monroe  desk  cal¬ 
culator  by  Robert  W.  Rockwell.  The  inverse  matrix  involved  contained  683  zeros. 

The  Gauss-Seidel  method  was  used  to  invert  the  matrix. 

Table  II  and  Figures  1  and  2  summarize  the  six-group  %n/  n  ratios  for  various 

values  of  k. 
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VII.  PROPAGATION  OF  PRECISION  INDEXES 


The  propagation  of  precision  indexes  was  accomplished  for  three  values  of  k  for 
Case  2  in  order  to  obtain  a  feeling  for  the  magnitude  of  the  uncertainty  in  the  variance/ 
mean  ratio.  Following  the  procedures  given  on  pages  205-208  of  reference  7  and  using 
the  constants  of  (74),  (75)  and  (76),  the  following  results  were  obtained  and  are  tabu¬ 
lated  below  in  Table  I  and  shown  on  Figure  2. 

Table  I 

Theoretical  Variance/Mean  Ratios  ±  Standard  Deviations 
of  the  Neutron  Population  for  Three  Values  of  k  for  Case  2 
(One  Delayed  Neutron  Group) 


k 

variance/mean  ratio  +  standard  deviation 

0.93 

1.62  ±  0.18 

0.97 

2.99  + 0.41 

0.995 

8.81 ±  1.38 

VIII.  NUMERICAL  RESULTS 

Table  II  and  Figures  1  and  2  tabulate  and  plot  the  theoretical  variance/mean 
ratios  of  the  neutron  population  in  the  subcritical  condition  of  the  AFRRI-TRIGA 
reactor. 
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Table  II 


* 


Theoretical  Variance /Mean  Ratios  of  the  Neutron  Population  in  the 
Subcritical  Condition  of  the  AFRRI-TRIGA  Reactor 


variance/mean  ratio 

Case  1* 

Case  2** 

Case  3*** 

0 

1.000 

1.000 

_ 

0.10 

0.951 

0.951 

— 

0.20 

0.905 

0.905 

— 

0.30 

0.863 

0.862 

— 

0.40 

0.827 

0.825 

— 

0.50 

0.800 

0.797 

— 

0.60 

0.790 

0.785 

— 

0.70 

0.813 

0.804 

— 

0.80 

0.920 

0.898 

— 

0.85 

1.056 

1.019 

1.037 

0.86 

1.098 

1.055 

1.074 

0.87 

1.  147 

1.098 

1.118 

0.88 

1.205 

1.148 

1.171 

0.89 

1.274 

1.208 

1.233 

0.90 

1.359 

1.281 

1.307 

0.91 

1.464 

1.370 

1.400 

0.92 

1.597 

1.480 

1.515 

0.93 

1.769 

1.622 

1.661 

0.94 

2.001 

1.807 

1.853 

0.95 

2.328 

2.061 

2.116 

0.96 

2.822 

2.426 

2.493 

0.97 

3.649 

2.995 

3.081 

0.98 

5.308 

3.999 

4. 120 

0.99 

10.297 

6.231 

6.429 

0.995 

20.285 

8.814 

9.106 

0.9995 

200.1 

14.324 

14.918 

0.99995 

1998 

15.295 

16.95 

0.999995 

19982 

15.400 

17.20 

Case  1  (no  delayed  neutron  groups)  was  computed  from  equation  (16). 

Case  2  (one  delayed  neutron  group)  was  computed  from  equation  (43). 

Case  3  (six  delayed  neutron  groups)  was  computed  from  equations  (44)  through  (71). 
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Figure  1.  Theoretical  variance/mean  ratio  of  the  neutron  population 
in  the  subcritical  condition  of  the  AFRRI-TRIGA  Reactor. 


Figure  2.  Theoretical  variance/mean  ratio  of  the  neutron  population 
in  the  subcritical  condition  of  the  AFRRI-TRIGA  Reactor. 
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